Is BFKL ruled out by the Tevatron gaps between jets data? 
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We have performed a detailed phenomenological investigation of the hard colour singlet exchange process 
observed at the Tevatron in events that have a large rapidity gap between outgoing jets. We include the effects of 
multiple interactions to obtain a prediction for the gap survival factor. Comparing the data on the fraction of gap 
events with the prediction from BFKL pomeron exchange we find agreement provided that a constant value of 
a a is used in the BFKL calculation. Moreover, the value of a s is in line with that extracted from measurements 
made at HERA. 
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1. Introduction 

Events with large rapidity gaps in the hadronic 
final state and a large momentum transfer across 
the gap, characterised by the presence of a hard 
jet on each side of the gap, have been observed 
in both pp collisions at the Tevatron [ [l], ||, || [| 
and in jp collisions at HERA [ ||, || . Such events 
are unexpected in standard Regge phenomenol- 
ogy since the cross section is predicted to fall as 
~ s _Q l*l, where a ~ 0.25 GeV -2 , whilst events 
with |i| > 1000 GeV 2 are routinely observed at 
the Tevatron. Clearly some other explanation 
must be sought. Uniquely in diffractive physics, 
high-t events are amenable to the use of pertur- 
bative QCD since the gap producing mechanism 
is squeezed to small distances [ f?J. Such calcu- 
lations have been carried out within the lead- 
ing logarithmic approximation of BFKL [ || by 
Mueller and Tang [ |[, and it is the aim of this 
talk to present comparisons of these calculations 
with the latest data from the Tevatron. The sit- 
uation is greatly complicated by the possibility 
that rapidity gaps formed by whatever process 
can be destroyed by multiple interactions between 
spectator partons in the colliding hadrons. De- 
tailed comparisons made and conclusions drawn 
from any dynamic model of high-t rapidity gap 
formation must therefore include a careful treat- 
ment of such physics. In this analysis, we use a 
model implemented in the PYTHIA Monte Carlo 
generator to simulate the effects of multi-parton 
interactions. 



2. D0 data versus the BFKL pomeron 

The analysis presented here was stimulated to 
some extent by the recent D0 measurements [ ^| 
of the fraction of dijet events containing a large 
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Figure 1. D0 data compared with a BFKL cal- 
culation. Plot from [ 01 . 



rapidity gap as a function of Et2, the Et of the 
second hardest jet, and the rapidity difference be- 
tween the two leading jets, Arj. The D0 results 
are shown in figure |l|. Jets are found using a cone 
algorithm [ [llj with cone radius 0.7 and the 
OVLIM parameter set to 0.5. The inclusive dijet 
sample is defined by the following cuts: 



• l^i I) l^l > L9, i.e. jets are forward or back- 
ward 

• V1V2 < 0, i.e. opposite side jets 

• E T2 > 15 GeV 

• Arj > 4, i.e. jets are far apart in rapidity. 
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The sub-sample of gap events is obtained by em- 
ploying the further cut that there be no particles 
emitted in the central region \r)\ < 1 with energy 
greater than 300 MeV. The BFKL curve is clearly 
ruled out by the data. The D0 BFKL curves are 
based on the calculation of Mueller and Tang im- 
plemented into the standard HERWIG 5.9 release 
[ |l2[ ^2|. In particular, the asymptotic cross- 
section of [ [| is used; in the limit y = Arj ^> 1, 



da(qq -> qq) 
dt 



(C F a s )' 



2tt 3 exp(2o7oy) 
~ (7a s C A ((3)yr 



(1) 



where cuq — w(0) = C y 4(41n2/7r)a s . The a A s in 
the pre-factor runs with —t according to the two- 
loop beta function, ujq = 0.3 and the a s in the 
denominator = 0.25. The falling of the BFKL 
curve with increasing Et2 is driven by the run- 
ning of the coupling in the pre-factor since the 
gap fraction goes like ~ a^/aj,. 
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3. Key issues 



Figure 2. Jet r\ profiles 



In this analysis, we choose somewhat different 
parameters. We also use the full Mueller Tang 
calculation without the asymptotic approxima- 
tion. This is also available in HERWIG 5.9 [ [|| 
and is available from the authors. We choose to 
fix a s — 0.17. To leading logarithmic accuracy 
a s is simply an unknown parameter. Higher or- 
der corrections will indeed cause the coupling to 
run, however it is not clear how this should be 
done in a consistent way. In this paper we restrict 
ourselves to the leading logarithmic approxima- 
tion and treat the coupling as a free parameter. 
Moreover, we are guided by recent HERA data on 
the double dissociation process [ [li| which can be 
described by the leading logarithmic BFKL for- 
malism with a s = 0.17. We also note that a fixed 
coupling constant was needed in order to explain 
the high-i data on pp elastic scattering via three 
gluon exchange [ 14 1. Furthermore, NLO correc- 



tions suggest a fixed value for the leading eigen- 
value of the BFKL equation, w(0), [ [l5| which in 
turn suggests the use of a fixed coupling in the 
LLA kernel. 

Underlying events and gap survival 

As mentioned above, it is critical in any esti- 
mate of gap formation rates to take into account 
the possibility that gaps can be destroyed by sec- 
ondary scatters, which may be perturbative or 
non-perturbative, between spectator partons in 
the colliding hadrons. Several models are avail- 
able [ |l6], [17, 18 1, but it would be fair to say that 
all are as yet in a early stage of development and 
are not tuned to pp data. We choose the model 



as implemented in PYTHIA 6.127 [0. Here the 
probability to have several parton-parton inter- 
actions in the same collision is modelled using 
perturbative QCD. The probability for additional 
interactions is not fixed but varies according to 
an impact-parameter picture, where central col- 
lisions are more likely to have multiple interac- 
tions. The partons in the proton are assumed 
to be distributed according to a double- Gaussian 
as described in [ [l9| [Tt]]. There are several pa- 
rameters in this model and we have used the de- 
fault setting for each.[] Our strategy is to gener- 
ate high-i photon exchange events (hard BFKL 
pomeron exchange has not been implemented in 
PYTHIA ) with and without multiple interactions, 
and take the percentage change in the number of 
rapidity gap events, defined as in the D0 analy- 
sis, as the gap survival factor. We find that gap 
survival in this model is to first order indepen- 
dent of E 3 ^ and A77, i.e. it can be treated as 
a multiplicative factor. The gap survival factor 
S does vary strongly with centre of mass energy, 
which is not unexpected since the number density 
of partons in the colliding hadrons, and there- 
fore the probability of having a secondary scat- 
ter, increases with energy. In summary, we find 
5(1800 GeV) = 22%, 5(630 GeV) = 35%. Full 
details can be found in [ g(J ■ 

A key point to notice is the interplay be- 
tween gap survival and underlying event : mul- 

1 Setting the switch MSTP(82)=4 in PYTHIA , with every- 
thing else default, will give the model as we have used 
it. 



2 



tiplc interactions also give rise to the so-called 
jet pedestal and underlying event effects. This 
means that the jets measured in hadron-hadron 
collisions cannot be compared directly to e.g. 
predictions from fixed order perturbation theory. 
In Figure [2] we show jet profiles obtained from 
PYTHIA with (mi4) and without (miO) multiple 
interactions (and with \5<j)\ < 0.7). The proton 
remnant is at (Jyy > 0. It is clear that multiple in- 
teractions introduce a jet pedestal of more than 
1 GeV of Et per unit rapidity. For comparison, 
also shown is the jet pedestal from HERWIG . We 
note that HERWIG predicts a greater amount of 
energy outside the jet cone than PYTHIA without 
multiple interactions. Again, a full discussion of 
these differences can be found in [ [2Cj . 

In the D0 jet measurements the excess Et 
from the underlying event is taken into account 
by correcting the jet Et using minimum bias 
data. In particular, the correction is determined 
by looking at the Et flow in regions away from 
the jets. The correction is made by subtracting 
approximately 1 GeV from the Et of each re- 
constructed jet [ In particular, in the gap 
fraction measurement, this subtraction is per- 
formed for all jets, including those in gap events. 
But, requiring a large rapidity gap also selects 
events without multiple interactions, where the 
jet pedestal is absent, or at least much smaller; 
multiple interactions destroy gaps, and therefore 
a gap event cannot have a multiple interaction. 
Since jet cross sections fall faster than 1/E T , such 
a correction can decrease the measured jet rate by 
up to 30% for 18 GeV jets. Our contention there- 
fore is that the jets in gap events should not be 
corrected for underlying event, and therefore the 
gap fraction should rise less steeply with Et than 
in figure [lj. 



4. Gap fractions 

Figures || and || show our results for the gap 
fractions as functions of An and Eti respectively. 
The stars are the HERWIG BFKL simulation with 
fixed a s — 0.17, with 1 GeV subtracted from each 
jet in order to simulate the D0 underlying event 
correction and the open circles are the D0 data. 
The gap fractions are constructed using a stan- 
dard PYTHIA QCD simulation without colour sin- 
glet exchange, and without multiple interactions. 
We have used both CTEQ2M and CTEQ3M par- 
ton distribution functions [ |24|, [2fj], and have 
found the differences to be small. Our philoso- 
phy is that the D0 data have been corrected for 
the effects of multiple interactions in non-singlet 
exchange events, and we should therefore gener- 
ate none, whereas we must undo the erroneous 
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Figure 3. Gap fraction as a function of At? com- 
pared to the D0 data 



* bfkl fixed <x s - 1 .0 GeV / pythia miO 
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Figure 4. Gap fraction as a function of Et2 com- 
pared to the D0 data 



correction to the colour singlet sample. The com- 
bination of fixing a s and correcting the gap events 
erroneously for multiple interactions produces the 
rise of the gap fraction at low jet Et- The solid 
circles show the gap fraction using a running a s 
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Figure 5. Gap fraction compared to the CDF 
data 



gap fractions at 630 GeV and 1800 GeV. We 
find that, even including gap survival effects, 
i?(630/1800) - 1 at the parton level. When 
hadronisation effects are taken into account, how- 
ever, we find that the ratio rises significantly to 
~ 3, with a strong dependence on Arj. D0 find 
i?(630/1800) = 3.4 ± 1.2 [ §, and CDF find 
i?(630/1800) = 2.4±0.9. In the D0 case the effect 
may be attributed to the different parton x ranges 
of the 630 GeV and 1800 GeV measurements (al- 
though we note that the CDF result is calculated 
at fixed x). The restriction x < 1 forces the gap 
and non-gap cross-sections to fall to zero at some 
maximum Ar/, Ar) ma]t . Now, the colour connec- 
tion that exists between the jets in the non-gap 
sample drags the jets closer together in rapid- 
ity. This has a small effect away from Ar7 max 
(since the Ar] spectrum is roughly flat) however 
as A77 — > A77 max it leads to a more rapid vanish- 
ing of the non-gap cross-section than occurs in the 
gap cross-section. This effect, combined with the 
fact that A77 max (630 GeV) < Ar? max (1800 GeV), 
leads to an enhancement of the measured 630 
GeV gap fraction at large Arj at the hadron level, 
and hence the larger value of i?(630/1800). 



in the BFKL sample. Even with the underlying 
event correction, this sample is unable to fit the 
data. The overall normalisation of the simulated 
gap fractions is multiplied by a factor of 0.6. That 
this is a reasonable thing to do can be appreci- 
ated once it is realised that our results have not 
been fitted to the data and that the overall nor- 
malisation is acutely sensitive to the magnitude 
of a s . Furthermore, the overall normalisation of 
the BFKL cross-section is uncertain since, within 
the leading logarithmic approximation, one does 
not know a priori the scale at which to evaluate 
the leading logarithms. Given these points, we 
conclude that the D0 data are in agreement with 
the leading order BFKL result. Figure || shows 
our result for the gap fraction as a function of 
Arj = 2rj* compared to the CDF data [ |j. Note 
that CDF do not attempt to correct their jets 
to include the effect of an underlying event. We 
therefore generate the PYTHIA non-singlet sample 
with multiple interactions (labelled mi4), and do 
not perform the 1 GeV / jet subtraction from the 
HERWIG BFKL sample. In this plot, our theory 
points are obtained using a renormalisation fac- 
tor of unity (compared to 0.6 in the D0 case). 
We then find reasonable agreement with the data 
except at the larger values of 77* where we are 
quite unable to explain a fall in the rj* distribu- 
tion. Recall however that D0 do not see a fall 
at large Ar]. Further clarification of the situation 
will require an increase in statistics. 

We have also computed the ratio of the 



5. Conclusions and future possibilities 

We have explicitly demonstrated that the Teva- 
tron data on the gaps-between-jets process at 
both 630 GeV and 1800 GeV are in broad agree- 
ment with the predictions obtained using the 
leading order BFKL formalism. However, we are 
not able to explain the behaviour of the CDF gap 
fraction at large A77. Agreement is obtained using 
the same fixed value of a s — 0.17 as was used to 
explain the recent HERA data on high-t double 
diffraction dissociation. 

Care must be taken in the interpretation of 
our findings, however. The BFKL formalism it- 
self suffers from being evaluated only to leading 
logarithmic accuracy. The uncertainties of the 
overall normalisation which follow will not be re- 
moved until an understanding of BFKL dynamics 
at non-zero t beyond the leading logarithmic ap- 
proximation in achieved. 

An understanding of the effects of underlying 
event and its impact on gap survival is crucial to 
the interpretation of the gaps between jets data, 
and indeed diffractive data as a whole. 

As pointed out in [ ^, |2(J , the gap fraction de- 
fined in terms of a region void of hadronic ac- 
tivity is not strictly infrared safe. A better ob- 
servable would be to define a gap to be a region 
that does not contain any jets with transverse 
momenta above some perturbatively large scale. 
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Work along these lines has also been performed 
in [ |6). 

One major disadvantage of the gaps between 
jets process arises from the need to measure both 
jets since this limits the reach in rapidity. In [ 
p7| , it was suggested to focus instead on the dou- 
ble dissociation sample (the gaps between jets 
events form a subsample of this generally much 
larger sample). By dropping the requirement to 
observe jets one not only gains in rapidity reach 
and statistics but also from the reduced systemat- 
ics associated with this more inclusive observable. 
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